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The bottom part of the neutron star crust is investigated using the Skyrme-Hartree-Fock 
approach with the Coulomb interaction treated beyond the Wigner-Seitz approximation. 
A variety of nuclear phases is found to coexist in this region. Their stability and relative 
energies are governed by the Coulomb, surface and shell energies. We have also found 
'"q ■ that a substantial contribution is coming from the spin-orbit interaction. 

The crust of a neutron star with a typical mass of 1.4 M , contains only 1.2% of the 
total stellar mass. Nevertheless its structure is important for understanding of several 
^ observational issues and thus provides a motivation for theoretical studies (see e.g. 
£ : § S i I I i % % i 0, and references therein). 

According to various theoretical models in the inner part of the crust, due to high 
density and pressure, the nuclei forming a crystalline lattice are immersed in a neutron 
gas. The electrons at these densities are ultrarelativistic particles and are assumed to 
form a uniform background. 

In the bottom part of the inner crust the interplay between the Coulomb and the surface 
energies leads to the appearance of various nuclear phases, characterized by different 
shapes. The early predictions, based on the liquid drop or Thomas-Fermi models, with 
quite restrictive symmetry conditions, lead to the considerations of five phases, formed in 
the region where the nucleon density varies from 0.03 to 0.1 fm~ 3 . These are: spherical 
nuclei, rods ("spaghetti" phase), slabs ("lasagna" phase), tubes and bubbles (see [ and 
references therein). 

These approaches however provided an oversimplified description of the system. First, 
due to the Wigner-Seitz approximation which has been used, there is no distinction be- 
tween phases resulting from different lattice geometries, second, the quantum effect as- 
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sociated with the shell correction energy has been neglected, and last but not least, the 
symmetry conditions, which admit only spherical, cylindrical and planar structures, are 
not sufficient to describe properly the variety of nuclear shapes present in the crust. 

Our approach, which is free from the aforementioned limitations, can shed some light on 
the complexity of the inner crust structure. The ground and excited (isomeric) states of 
neutron-proton-electron (npe) matter have been determined by the Skyrme-Hartree-Fock 
method. We have solved the HF equations by discretization in coordinate space, within a 
cubic box and with periodic boundary conditions imposed on the nucleon wave functions 
T2|j . The details of this approach have been described in Ref. [ 



In the HF approximation, the total energy of a nuclear system is, in general, a functional 
of three local densities: total nucleon density p(r) = Pp(r) + p„(r), total kinetic density 
r(r) = Tp(r) + r n (r) and total spin-orbit density J(r) = J p (r) + J n (r). Each of the densities 
is in turn expressed through the single-particle wave functions 4>f (r, a): 

P 9 W=E^I^(r,a)| 2 , (1) 



nl\ V0?(r,a)| 2 , (2) 



J,(r) = -i £ <0f(r, orjVtf (r, a') x (a | a \ a') (q = n,p), (3) 



where \a (a = ±1) is the spin projection on the z-axis, and n\ is the occupation factor, 
which is 1 below the Fermi level and zero otherwise. 



E 



The Skyrme energy functional can be written as 
^ 2 



—t + B x p + B 2 (p p + p n ) + B 3 pr + B 4 (p p T p + p n T n ) + 



2m 

+ B 5 pV 2 p + B 6 (p p V 2 p p + p n V 2 Pn ) + 5 7 p 2+7 + 

+ Bsp^ipl + pi) + Bg( P V ■ J + p p V • J p + PnV ■ J r 



d s r. (4) 



The coefficients Bi depend on the Skyrme force parameters [ |12| . It is useful to distinguish 
the following contributions to the Skyrme energy functional: 



Egksurf 



E x 



J (B 5 pV 2 p + B 6 ( Pp V 2 Pp + p n V 2 p n )) d 3 r, (5) 
J (5 9 (pV • J + ppV ■ Jp + p„V • J n )) d 3 r, (6) 



where E s k SUT f contains the contribution to the surface energy and E so is the spin-orbit 
term. 

For a correct description of npe matter one has to add to the Skyrme energy functional 
the Coulomb term and the kinetic energy of electrons: 



E — E s k + Ecoui + E e i, 



(7) 
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where the direct Coulomb energy has been calculated by solving the Poisson equation 
with periodic boundary conditions. The Coulomb exchange term for both protons and 
electrons has been treated in the Slater approximation. 

Typical results for two values of the average total density p av = — / (p p + p n )d 3 r, and 

V Jv 

the proton-to-nucleon ratio (Z/A) ensuring the (3 stability condition, are shown in Fig. 
1 . The total energy densities have been plotted with respect to the value of the spherical 
system Q p = 0. 
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Figure 1. The results of the Skyrme-Hartree-Fock calculations (using SLy4 force) pre- 
sented as a function of the axial proton quadrupole moment: Q20 — Qp- The calculations 
have been performed in a box of lengths: d = 26/m (left panel) and d = 20.8/m (right 
panel). The upper subfigures show the total energy density and the lower subfigures the 
contribution to the total energy coming from various terms in the density functional (see 
text for details). The phase A consists of spherical nuclei in the sec lattice, the phases 
B and C denote the sec geometries differing by the nuclear deformations. The D phase 
corresponds to the bec lattice. 

The energy density has several local minima, corresponding to various nuclear shape 
and lattice geometries. It was shown in Ref. [ [K| that the relative position of various 
nuclear configurations is a sensitive function of the density. It is a consequence of the 
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shell effects associated with unbound neutrons [ fU|]. Hence different nuclear phases may 
coexist in the same density range. There is an obvious larger variety of nuclear shapes 
than it was predicted in the previous approaches ([ 0, and references therein). 

In the lower subfigures of Fig. 1, the contributions from the Coulomb energy, the E s k sur f 
term and the spin-orbit term, have been shown. One can make the following observations: 

• The Coulomb energy decreases as a function of Q p for small values of quadrupole 
moment and then increases slightly for large Q p values. This is related to the 
fact that large Q p values correspond to the configuration where the nuclei from 
neighboring cells join each other. Hence the configuration favored by the Coulomb 
interaction consists of deformed nuclei, but not necessarily in the shape of rods or 
slabs. 

• The part of the surface energy coming from the E s k SU rf term follow approximately 
the behavior of the term E s f~ — E so . It suggests that the surface energy dominates 
the behavior of E.h function of the deformation. 

• The amplitude of variation of the spin-orbit term E so is of the same order as the 
E s k — E so . The spin-orbit energy is a sensitive function of Q p and clearly have a 
direct influence on the appearance and depths of local minima of the total energy 
curve. One has to remember that this term gives the substantial contribution to 
the shell energy. 

Finally one has to mention that the pairing correlation were not included in the above 
calculations. The pairing would smooth the total energy curve and thus its role in the 
inner edge of the neutron star crust needs a thorough examination. 
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